Optical spectroscopic study of the interplay of spin and charge in a'-NaV20E 
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We investigate the temperature dependent optical properties of a'-NaV205, in the energy range 
4ineV-4eV. The symmetry of the system is discussed on the basis of infrared phonon spectra. By 
analyzing the optically allowed phonons at temperatures below and above the phase transition, we 
conclude that a second-order change to a larger unit cell takes place below 34 K, with a fluctuation 
regime extending over a broad temperature range. In the high temperature undistorted phase, we 
find good agreement with the recently proposed centrosymmetric space group Pmmn. On the other 
hand, the detailed analysis of the electronic excitations detected in the optical conductivity, provides 
direct evidence for a charge disproportionated electronic ground-state, at least on a locale scale: A 
consistent interpretation of both structural and optical conductivity data requires an asymmetrical 
charge distribution on each rung, without any long range order. We show that, because of the locally 
broken symmetry, spin-flip excitations carry a finite electric dipole moment, which is responsible 
for the detection of direct two-magnon optical absorption processes for E\\a. The charged-magnon 
model, developed to interpret the optical conductivity of a'-NaV205, is described in detail, and its 
relevance to other strongly correlated electron systems, where the interplay of spin and charge plays 
a crucial role in determining the low energy electrodynamics, is discussed. 



I. INTRODUCTION 

After many years of intensive experimental and theo- 
retical work on CuGeOa another inorganic compound, 
a'-NaV2 05, has been attracting the attention of the 

■ scientific community working on low-dimensional spin 
[ systems, in general, and on the spin-Peierls (SP) phe- 
• nomenon, in particular. In fact, in 1996 the isotropic ac- 
tivated behavior of the magnetic susceptibiliiv has been 

. observed at low temperatures on a'-NaV205E Moreover, 
' superlattice reflections, with a lattice modulation vec- 
tor /c = (7r/a,|jr/&, 7r/2c), were found in an x-ray scattering 

■ experimentjD and a spin gap A=9.8 meV was observed, 
] at the reciprocal lattice point (27r/a, ir/b, 0), with inelastic 

■ neutron scattering.0 The SP picture was then proposed 
] to explain the low temperature properties of this cont 
. pound, with a SP transition temperatures Tsp=34 KB 
' In this context, a'-NaV205 seemed to be a particularly 

interesting material, as far as the full understanding of 
the SP phenomenon is concerned, because on the basis 
, of the structural analysis reported by Carpy and Galy 
1 in 1975,_suggesting the noncentrosymmetric space group 
P2imnp it was assumed to be described by the one- 
dimensional (ID) spin-1/2 Heisenberg (HB) model bet- 
ter than CuGeOs. In fact, for the latter compound it has 
been established both experimentallyQ and theoreticallyQ 
that the 2D character of the system cannot be neglected. 
On the other hand, on the basis of experimental re- 
sults later obtained, the interpretation of the phase tran- 
sition is still controversial. First of all, the reduction 
of Tsp upon increasing the intensity of an externally 
applied magnetic field, expected for a true SP system. 



has not been observed. Q'lI For instance, Biichner et al.JB 
performing low temperature magnetization and specific 
heat measurements in magnetic field if < 14T, obtained 
ATsp(H) = rsp(0)-rsp(14T)«0.15K, i.e., almost a fac- 
tor of 7 smaller that-pYhat expected on the basis of Cross 
and Fisher theory.™ Furthermore, they claimed that 
the entropy reduction experimentally observed across the 
phase transition is considerably larger than the theo- 
retical expectation for a ID antiferromagnetic (AF) HB 
chain with only spin degrees of freedom and J = 48 meV 
(i.e., value obtained for a.'-NaV205 from magnetic sus- 
ceptibility measurementscl) . In order to explain these 
findings, degrees of freedom additional to those of the 
spin system should be taken into account, contrary to the 
interpretation of a magnetically driven phase transition^ 
This already quite puzzling picture became even more 
complicated when the crystal structure and the sym- 
metry of the electronic configuration of a'-NaV205, in 
the high temperature phase, were investigated in de- 
tail. In factj—oa, the basis of new x-ray diffraction 
measuremcntslljil3 it was found that the symmetry of 
this compound at room temperature is better described 
by the centrosymmetric space group Pmmn than by the 
originally proposed P2imnB On the other hand, the in- 
tensities and polarization dependence of the electronic 
excitations detected in optical conductivity spectra gave 
a direcjt|-evidence for a broken-parity electronic ground 
state,Ej'L3 in apparently better agreement with the non- 
centrosymmetric space group P2imn. The solution of 
this controversy and the assessment of the symmetry is- 
sue are of fundamental importance for the understanding 
of the electronic and magnetic properties of the system 



and, ultimately, of the phase transition. 

In this paper we present a detail investigation of the 
interplay of spin and charge in a'-NaV205, on the basis 
of optical reflectivity and conductivity data in the energy 
range 4 nieV-4 eV. We will start from the discussion (fttte 
newly proposed high temperature crystal structurejljilj 
and its implications on the physics of this compound (sec- 
tion |l|). A group theoretical analysis for the two differ- 
ent space groups will be reported and later compared to 
the experimentally observed far-infrared phonon spectra. 
In section 01, we will present a theoretical model which 
we will use to interpret qualitatively and quantitatively 
the optical conductivity spectra of a'-NaV205. In sec- 
tion [V, we will analyze the data obtained with optical 



spectroscopy at different temperatures above and below 
Tc = 34 K (throughout the paper, we will refer to the 
transition temperature as Tc because the interpretation 
of the phase transition is still controversial). In particu- 
lar, we will concentrate on: 

(i) Analysis of the temperature dependent phonon 
spectra, in relation to the x-ray diffraction data, in order 
to: First, learn more about the symmetry of the mate- 
rial in both high and low temperature phase. Second, 
detect signatures of the lattice distortion and study the 
character of the phase transition. 

(ii) Detailed study of the very peculiar electronic and 
magnetic excitation spectra. In particular, we will show 
that we could detect a low frequency continuum of ex- 
citations which, on the basis of the model developed in 
section [II, we ascribe to 'charged bi-magnons', i.e., di- 
rect two-magnon optical absorption processes. 

Finally in section M, we will discuss the relevance of our 
findings to the understanding of the nature of the phase 
transition in a'-NaV205. In particular we will try to as- 
sess whether the picture of a charge ordering transition, 
accompanied by the opening of a magnetic gap, is a valid 
alternative to the originally proposed SP description. In 
fact, this_alt^native interpretation has been recently put 
forward,Ej~E3 on the basis of temperature dependent nu- 
clear magnetic resonance (NMR) data.Ej 




FIG. 1. High temperature crystal structure of a'-NaV205 
Pyramids around V ions and rows of Na ions are shown. 



grouped into sets of two, forming a ladder, with the rungs 
oriented along the a axis. The rungs are formed by two V 
ions, one on each leg of the ladder, bridged by an O ion. 
The V-0 distances along the rungs are shorter than along 
the legs, implying a stronger bonding along the rung. In 
the a-b plane the ladders are shifted half a period along 
the b axis relative to their neighbors. The noncentrosym- 
metric space group P2imn, allows for two inequivalent V 
positions in the asymmetric unit. These sites were inter- 
preted as different valence states, V^"*" and V^"*", repre- 
sented by dark and light gray balls, respectively, in Fig. 
0. In this structure it is possible to identify well-distinct 
ID magnetic V''+ (S=l/2) and non-magnetic V^+ (S=0) 
chains running along the b axis of the crystal, and al- 
ternating each other along the a axis. This configura- 
tion would he responsible for the ID character (Bonnee, 
Fisher-typa23) of the high temperature susceptibilityE^I 
and for the SP transitioH, possibly involving dimerization 
within the Y'^'^ chains.lil Moreover, due to the details of 
the crystal structure, only one of the thcee t2g orbitals 
of V''"'", i.e., the dxy orbital, is occupiedJla As a result, 
the insulating character of a'-NaV205, ■wjhJA has been 
observed in DC resistivity measurements,c3a would be 
an obvious consequence of having a ID 1/2-filled Mott- 
Hubbard system. 



II. THE SYMMETRY PROBLEM IN a'-NAVzOs 

The interpretation of the phase transition at 34 K as a 
SP transitiorcl is based on the noncentrosymmetric space 
group P2imn originally proposed by Carpy and Galy for 
the high temperature phase of Q;'-NaV205.l3 Following 
their crystallographic analysis, the structure can be con- 
structed from double-rows (parallel to the b axis) of edge 
sharing pyramids, one facing up and the other down, with 
respect to the a-b plane (see Fig. |^). These double rows 
are connected by sharing corners, yielding a planar ma- 
terial. The planes are stacked along the c axis, with the 
Na in the channels of the pyramids. In the a-b plane (see 
Fig. 0) , we can then identify linear chains of alternating 
V and O ions, oriented along the b axis. These chains are 




FIG. 2. Two-leg ladder structure formed by O and V ions in 
the a-b plane of Q'-NaV205. Black dots are O ions; dark and 
light gray balls are V*''' and V^"'' ions, respectively, arranged 
in ID chains in the space group P2\mn. 



Howeveij, |_recent redeterminations of the crystal 
structurell3il3 showed that the symmetry of a'-NaV205 
above T^ is better described by the centrosymmetric 
space group Pmmn: The suggested structure is, very sim- 
ilar to the one proposed by Carpy and Galyu (a small 
difference is the somewhat more symmetric eight-fold co- 
ordination of NallJ). On the other hand, one major dif- 
ference was found, namely, the presence of a center of 
inversion in the unit cell. The evidence is the very low 
i?(-F)-value of 0.015, and the fact that no lower value of 
R{F) can be obtained when the x-ray diflfraction data 
are refined omitting the inversion center.c3 In particu- 
lar, this result suggests the existence of only one kind 
of V per unit cell, with an average valence for the V 
ions of -1-4.5. This finding has important implications for 
the understanding of the electronic and magnetic prop- 
erties of Q;'-NaV205. In fact, we have now to think in 
the framework of a 1/4-filled two-leg ladder system: It 
is still possible to recover an insulating ground state as- 
suming that the d electron, supplied by the two V ions 
forming a rung of the ladder, is not attached to a partic- 
ular V site but is shared, in a V-O-V molecular bonding 
orbital along the rungJlj In this way, we are effectively 
back to a ID spin system (a two-leg ladder with one spin 
per rung in a symmetrical position). However, in case of 
a SP phase transition, singlets would have to be formed 
by electrons laying in molecular-like orbitals and a rather 
complicated distortion pattern would have to take place: 
Not simply the dimerization within a real ID chain (as, 
e.g., in CuGeOs), but a deformation of the plaquettes, 
formed by four V ions, within a ladder. 

Regarding the controversy about the x-ray structure 
determination for a'-NaV205 above Tc, it is important 
to stress that, even though it is now well established that 
the 'most probable' space group is the centrosymmetric 
Pmmn (on the basis of the statistics of the refinement), 
deviations from centrosymmetcy, up to 0.029 A are un- 
likely but cannot be excluded.ll3 Moreover, one has to 
note that x-ray diffraction measurements are sensitive to 
the charge distribution of core electrons and not of va- 
lence electrons. Therefore, if a local breaking of symme- 
try in the distribution of the V d electrons were present on 
a local scale without the long range order required by the 
noncentrosymmetric space group P2imn (as suggested 
by th e optical conductivity data we will present in section 
IV D|) , it would not be detectable in an x-ray diffraction 
experiment: X-rays would just see the 'space-average' 
given by the centrosymmetric space group Pm,m,n. 



A. Group Theoretical Analysis 



of modes has to be expected because each unit cell con- 
tains two formula units, corresponding to 48 degrees of 
freedom, which will be reflected in 48 phonon branches. 
In the space group Pmmn the site groups in the unit 
cell are: C2y for the 2 Na and the 2 0(1) atoms, and 
Cf^ for the 4 V, the 4 0(2), and the 4 0(3). Follow- 
ing the pjiclcar site group analysis method extended to 
crystals,Eil the contribution of each occupied site to the 
total irreducible representation of the crystal is: 

rNa+0(l) —'^[-Ag + Blu + B2g + B2u+ B3g + B^y] , 

rv+o(2)+o(3) =3[2Ag + Au + Big + 2Biu 

+2B2g + B2u + B3g + 2B3u] . 

Subtracting the silent modes (SA^) and the acoustic 
modes {Biu + B2u + B3u), the irreducible representation 
of the optical vibrations, for the centrosymmetric space 
group Pmmn, is: 

r^ 8Ag{aa,bb,cc) + 3Big{ab) + 8B2g {ac} + bBsg (be) 
+7BiuiE\\c)+AB2u{E\\b) + 7B3uiE\\a) , (1) 

corresponding to 24 Raman {Ag,Big,B2g,B3g) and 18 in- 
frared {Biu,B2u,Bzu) active modes. 

In the analysis for the noncentrosymmetric space group 
P2imn, we have to consider only the site group C^^ for 
aU the atoms in the unit ceU: 2Na, 2V(1), 2V(2), 2 0(1), 
2 0(2), 2 0(3), 2 0(4), and 2 0(5). Therefore, the total 
irreducible representation is: 

8 X Tcy'^8x[2Ai + A2 + Bi + 2B2]. 

Once again, subtracting the acoustic modes Ai + Bi + 
B2 (there is no silent mode in this particular case), the 
irreducible representation of the optical vibrations, for 
the noncentrosymmetric space group P2imn, is: 



r'^15Ai{aa, bb, cc; £;||a) -1-8^2(60) 
+7Bi{ab- E\\b) + 15B2{ac- E\\c) , 



(2) 



corresponding to 45 Raman {Ai,A2,Bi,B2) and 37 in- 
frared {Ai,Bi,B2) active modes. 

As a final result, we see that the number of optical vi- 
brations, expected on the basis of group theory, is very 
different for the two space groups. Moreover, whereas in 
the case of Pmmn the phonons are exclusively Raman 
or infrared active, for P2imn, because of the lack of in- 
version symmetry, there is no more distinction between 
gerade and ungerade and certain modes are, in principle 
(group theory does not say anything about intensities), 
detectable with both techniques. 



An alternative way to assess the symmetry issue is to 
investigate the Raman and infrared phonon spectra, com- 
paring the number of experimentally observed modes to 
the number expected for the two different space groups 
on the basis of group-theory. Obviously, a large number 



III. CHARGED-MAGNON MODEL 

In the context of the extensive work done, in the 
last few years, on 2D S=l/2 quantum AF and their 
ID analogs, it has been shown both experimentally and 



theoreticallyETEJ that also optical spectroscopy, besides 
elective techniques like Raman and neutron scattering, 
can be a very useful probe in studying the spin dynamics 
in this systems. For example, in the 2D-system La2Cu04 
(ID-system Sr2Cu03) where two-magnon (two-spinon) 
excitations are in principle not optically active because 
of the presence of a center of inversion that inhibits the 
formation of any finite dipolc moment, phonon- assisted 
magne tic , e xcitations were detected in the mid-infrared 
region.EaG3 These magneto-elastic absorption processes 
are optically allowed because the phonon there involved 
is effectively lowering the symmetry of the system. 

Let us now concentrate on a system where a breaking 
of symmetry is present because of charge ordering. We 
show that, in this case, magnetic excitations are expected 
to be directly optically active and detectable in an opti- 
cal experiment as charged bi-magnons: Double spin-flip 
excitations carrying a finite dipole moment. We will do 
that investigating the spin-photon interaction in a single 
two-leg ladder with a charge disproportionated ground 
state. On the basis of this model we will later interpret 
the optical conductivity spectra of a'-NaV2 05. 



A. Model Hamiltonian 

In this section we discuss the single two-leg ladder de- 
picted in Fig. g, where t± and ty are the hopping pa- 
rameters for the rungs and the legs, respectively, d± is 
the length of a rung, and -A/2 and -I- A/2 are the on-site 
energies for the left and the right leg of the ladder, re- 
spectively. We work at quarter filling, i.e., one electron 
per rung. The total Hamiltonian Ht of the system is: 



Ht^ Ho + H±+ H» 



(3) 



Ho=uY^ {njR^rijRl + UjL^rijLi} + A ^ n^c (4) 
i i 

i7x=ii5]{4,i?,,. + i?.c.} (5) 



J,cr 



H\\ = h 



E {4.^^+1^- + lLl,+i.. + H.c] , (6) 



J.cr 



where L^ ^ (i?J ^) creates an electron with spin a on the 
left-hand (right-hand) site of the jth rung, U is the on- 
site Hubbard repulsion, njj^^^rijjii {njL^njLi) is counting 
the double occupancies on the right-hand (left-hand) site, 
njc = {njR — "'jl)/2 is the charge displacement operator, 
and A is the potential energy difference between the two 
sites. Considering only Hq at quarter filling one imme- 
diately realizes that A is not only breaking the left-right 
symmetry of the system (for the symmetric ladder A = 0) , 
but it is also introducing a long range charge ordering 
with one electron per site on the left leg of the ladder. 
Working in the limit t±_ ^i|| we can consider each rung 



as an independent polar molecule (with one electron) de- 
scribed by the Hamiltonian Hjq+Hj±. The two solutions 
are lob-sided bonding and antibonding wave functions 

\L) =u\L) + v\R), and \R) ^ u\R) - v\L) , with: 



u 


1 /i ^ 
" ^/2 V ' EcT ' 


ler 


e EcT = V^' + 4ii 



V2 



A 



E, 



CT 



(7) 



is the splitting between these 
two eigenstatcs. The excitation of the electron from the 
bonding to the antibonding state is an optically active 
transition with a degree of charge transfer (CT) from the 
left to the right site which is the larger the bigger is A. It 
is possible to calculate the integral of the real part of the 
optical conductivity CTi(ct') for this excitation, quantity 
that can be very useful in the analysis of optical spectra. 
We start from the general equation:Eil 



ai{Lu)duj 



n#9 



\{n\Y,^M\ (8) 



where V is the volume, qe is the electron charge and 
i is the site index. As for a single rung (i?|xi|L) = 
—d±t±/EcT, we obtain the following expression, which 
is exact for one electron on two coupled tight-binding 
orbitals: 



ai{LL))duj 



T:qlNd\t\h- 



2 77-1 



'-E. 



CT 



CT ' 



(9) 



where N is the volume density of the rungs. Comparing 
Eq. ^ and the expression for Ect with the area and the 
energy position of the CT peak observed in the experi- 
ment, we can extract both t± and A. 

In the ground state of the Hamiltonian Ho+H±, each 
electron resides in a \L) orbital, with some admixture 
of \R). Let us now introduce the coupling of the rungs 
along the legs, considering _ff|| and a small fragment of 
the ladder with only two rungs. We have now to take into 
account the spin degrees of freedom: If the two spins are 
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FIG. 3. Sketch of the single two-leg ladder. Circles repre- 
sent the ionic sites having on-site energies -A/2 and -l-A/2 
on the left and on the right leg of the ladder, respectively. 
t_L and i|| are the hopping parameters for the rungs and the 
legs, respectively, d± is the length of a rung, and j is the rung 
index. Also indicated is the relation, assumed throughout the 
discussion, between the different energy scales of the model. 



parallel, the inclusion of H\\ has no effect, due to the 
Pauli-principle. If they form a S' = state, the ground 
state can gain some kinetic energy along the legs by mix- 
ing in states where two electrons reside on the same rung. 
We start from the singlet ground state for the two rungs: 



\L1L2 



- -7^ ||-^iT-^2i) - \LiiL2-\)j 



(10) 



with eigenvalue i?L 1L2 = ^Ect, and we study the ef- 
fect of 7?|| using a perturbation expansion in t\\/ Ecti to 
the second order. However, H\\\LiL2) does not coincide 
with any of the singlet eigenstates of the system consist- 
ing of two electrons on one rung with the Hamiltonian 
Hjo+Hj±. In order to obtain the exchange coupling con- 
stant Jii between two spins on neighboring rungs, we have 
first to calculate the set of spin-singlet eigenstates for the 
two electrons on one rung, and then to evaluate the frac- 
tional parentages of _ff|| \L1L2) with respect to that set. 
The Hilbert subspace of two-electron spin-singlet states 
is spanned by the vectors: 

\LjL,)^l1^l]J0) (11) 

\L,R,)^±{lIrI-LIrI}\0) (12) 

|i?,i?,)=4^4jO). (13) 

On this basis the Hamiltonian of the jth rung is: 



Hjo + Hj± 



U - A V2t± 

V2t±_ V2t±_ 

V2t± U + A 



(14) 



A significant simplification of the problem occurs in the 
large U limit which, moreover, elucidates the physics of 
exchange processes between a pair of dimers in an ele- 
gant and simple way. We therefore take the limit C/— >cxd. 
The solutions for A = arer3 The (anti)-symmetric lin- 
ear combinations of \LjLj) and \RjRj) with energy U, 



and 



LjRj) with energy Elr-- 



-4tyC/^0. The triplet 



states have energy 0. For a general value of A, still obey- 
ing A ^ C/, the only relevant eigenstate for the calculation 
of the exchange processes between neighboring dimers is 
\LjRj), with energy Elu — O. The other two states with 
energy of order U are projected out in this limit. 

We can now go back to the problem of two rungs with 
one electron per rung, and consider the hopping along 
the legs treating 7f || as a small perturbation with respect 
to Ho+H± . We proceed by calculating the corrections to 
\L1L2) by allowing a finite hopping parameter t\\ between 
the rungs. Using lowest order perturbation theory, the 
ground state energy of a spin-singlet is: 



Eg,S=0 — —EcT — (-^1-^21^1 



1 



Hn+Hi+E, 



-Hn\LiL2) 



CT 



(15) 



We proceed by calculating the coefficients of fractional 
parentage by projecting H\\\LiL2) on the two-electron 
single eigenstates of a single rung: 



C^^ j+1 i^i 



{ 



va^-Q 



C^J) J-1 




FIG. 4. Pictorial description of the electric dipole moment 
associated with a single spin-fiip in the asymmetrical two-leg 
ladder. If spins are antiparallel (left-hand side), the charge 
of the j th electron is asymmetrically distributed not only 
over the j th rung but also on the nn rungs (with opposite 
asymmetry), because of the virtual hopping. If we flip its 
spin (right-hand side) , the j th electron is conflned to the j th 
rung, resulting in a reduced charge density on the right leg 
of the ladder. Therefore, there is a net dipole displacement 
(along the rung direction) of the antiparallel spin configura- 
tion compared to the parallel one. 

TJjilZiZz) = V2t|| Y^ {u^\LjLj)+v^\RjRj) 

+V2uv\LjRj)}. (16) 

Calculating the second order correction to £'liL2 for the 
S—Q state, in the limit where U ^00, and realizing that 
the correction for the 5*^1 spin configuration is zero, we 
obtain the exchange coupling constant Jii : 



Jii 



StUl 



[A2-h4i' 



2l3/2 



(17) 



B. Interaction Hamiltonians and Effective Charges 

Let us now consider three rungs of the ladder (see Fig. 
Eh to understand the role of spin-flip excitations in this 
system. Because of the energy gain J\\ , in the ground 
state we have antiparallel alignment of the spins (Fig. 
H, left-hand side). The jth electron is described by the 
lob-sided bonding wave function \Lj) and, because of 
the virtual hopping to the nearest-neighbor (nn) rungs, 
by states where two electrons reside on the same rung. 
Working in the limit f/ — > 00 these states have one elec- 
tron in the \L) and one in the \R) state on the same rung. 
Therefore, the charge of the jth electron is asymmetri- 
cally distributed not only over the jth rung (dark grey 
area in Fig. H, left) but also on the nn rungs (light grey 
area in Fig.^, left). However, in the latter case most 
of the charge density is localized on the right leg of the 
ladder. If we flip the spin of the jth electron (Fig. ^, 
right-hand side) no virtual hopping is possible any more 
to the nn rungs because of the Pauli principle. The elec- 
tron charge distribution is now determined only by the 
lob- sided wave function \Lj) (dark grey area in Fig. ^, 



right). As a result, there is a net dipole displacement of 
the antiparallel spin configuration compared to the paral- 
lel one: The spin-flip excitations carry a finite dipole mo- 
ment parallel to the rung direction (a axis). This dipole 
moment can couple to the electric field of the incident 
radiation in an optical experiment, and result in what 
we refer to as charged bi-magnon excitations. One may 
remark, at this point, that only double spin fiips can be 
probed in an optical experiment because of spin conser- 
vation. However, what we just discuss for a single spin 
flip can be extended to the situation where two spin flips, 
on two different rungs, are considered. 

The coupling to light with E\\d (rung direction) can 
now be included using the dipole approximation. The 
only effect is to change the potential energy of the \R) 
states relative to the \L) states. In other words, we have 
to replace A with lS.-\-qi,d_i_Ea^ where Ea is the component 
of the electric field along the rung. The coupling to the 
CT transition on the rungs is given by the Hamiltonian: 



HcT = qed±Eanc = qed±Ea ^ rijc 



(18) 



The Hamiltonian for the spin-photon coupling can be ob- 
tained taking the Taylor expansion of Jii with respect to 
Ea, and retaining the term which is linear in the electric 
field. Noticing that dEa = dA/qed±, we have: 



Hs — <lmd_LEahs — qmd^Ea ^ Sj ■ S- 



j 



'J+1 1 



(19) 



dJii 



dJ, 



where 9m = j- g^ — 9e -g^ is the effective charge involved 

in a double spin-flip transition, and hs = J2n ^j'^j+i- ^^ 
the limit U- ' - ■ 



>oo, the effective charge reduces to: 

3J||A 



= A2+4ti 



(20) 



where J\\ is given by Eq. n^. One has to note that for 
a symmetrical ladder, where A = 0, the effective charge 
qm vanishes, and the charged magnon effect disappears: 
The peculiar behavior of the spin flips we just described 
is an exclusive consequence of the broken left-right sym- 
metry of the ground state. On the other hand, the 
bonding-antibonding transition would still be optically 
active. However, the bonding-antibonding energy split- 
ting would be determined only by ij^, i.e., E'Qj, = 2t±. At 
this point we like to speculate on the role played by A 
in determining not only the optical activity of the dou- 
ble spin-flip excitations but also the insulating (metal- 
lic) nature of the system along the leg direction. For 
A = and for strong on-site Coulomb repulsion U, the 
energy cost to transfer one electron from one rung to the 
neighboring one is Eqj, = 2t±. Therefore, the system is 
insulating for 4ij| < E^j, = 2t± and metallic in the op- 
posite case. If we now start from a metallic situation 
4i|| >E'(jrp — 2t±, switching on A we can realize the con- 
dition 4i|| < EcT = -\/AM-4tY which would result in a 
metal-insulator transition and, at the same time, in op- 
tical activity for the charged bi-magnons. 



C. Spectral Weights 

In this section we will attempt to estimate what would 
be, in an optical experiment, the relative spectral weights 
for the CT and the charged bi-magnon excitations of the 
single two-leg ladder. In particular, we want to show that 
these quantities can be expressed in terms of properly 
defined correlation functions. Considering the integral 
of the optical conductivity, Eq. ^, and the interaction 
Hamiltonians Hct and Hs, we can write: 
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Let us first concentrate on the CT excitation which is a 



simpler problem because, as we discussed in section [II A 



on each rung we are dealing with one electron on two 
coupled tight-binding orbitals \L) and \R). Therefore, 
only one energy term {En— Eg) = Eqt has to be considered 
in the previous equation. We can then write: 
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= EcT {{g\nl\g) - \{g\nc\g)\^} = 
^ EcT {(nl) - (nc)'} , 



(23) 



where we took into account that the expectation value of 
nc over the ground state is a real number. From Eq.|2l| 
andP3, we obtain: 



where we have introduced the correlation function: 
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For the ground state of A'^ independent rungs per unit 
volume \g) — Y\j\Lj), we calculate gc — N[N ~ [N ~ 
l)'iu'^v'^~N{v'^~'V?)'^]. Because {v'^-u^)'^ ^l- 'iv?v'^ , it 
follows that gc = N{2uv)'^ . Therefore, gc — ^ if the two 
sites on a rung are completely independent: u= 1 and 
u = 0, for fj^=0 and A^O. On the other hand, gc has its 
maximum value gc = N when the rungs form homo-polar 
molecules: u — v— 1/a/2, for t±^0 and A = 0. Noticing 
that l/(2w)2 = 1-1- (A/2i_L)2 one can see that Eq. || and 
Eq. II correspond exactly to the same result. 

In performing a similar quantitative analysis for the 
spin-flip excitations in the single two-leg ladder, we have 
to consider the elementary excitations of the ID S=l/2 
HB AF, system which does not have long range magnetic 



order at any temperature due to its ID ckaracter. It has 
been shown by Faddeev and TakhtajarO that the true 
elementary excitations of the ID HB AF are doublets of 
S=l/2 spinon excitations. Each spinon is described by 
the dispersion relation: 



e (k) — — ■ Jsinfc . 



< fc < TT. 



(26) 



The two-spinon continuum, with total spin S=l or S=0, 
is defined by e'(g) = e'(fci) + e'(fc2), where q = ki+k2, 
and fci and k2 are the momenta of the spinous. The 
lower boundary is found for ki = and k2 = q (or the 
other way around): e[{q) — ^Jsinq. The upper one for 
ki=k2=q/2: £2(9) = 7rJsin(q/2). 

In the present case, we are dealing with double spin 
flips generated by the interaction Hamiltonian Hs- What 
is then relevant is the four-spinon continuum defined by 
e'(g) =J2i c'(^j)i where q = J2i ^i^ ^^^ ^^^ index i labels 
the four spinous. Of this excitation spectrum we would 
actually probe in an optical experiment only the q = (or 
equivalently q = 2tt) states, with total spin S=0. These 
states form a continuum at g = extending from e'l (0) — 
to e2(0) = 27rJ. In fact, because Eq. |2^ has its maximum 
^max = 7rJ/2 for k = 7r/2, four spinous with the same 
k — tt/2 correspond to the q = 2tt {q = 0) excited state 
with the highest possible energy: e2(0) = 4e^(ij. = 27rJ. 

If we now go back to Eq. p2l we see that the evaluation 
of the integrated optical conductivity for the spin excita- 
tions is problematic because it requires a detailed analy- 
sis of the matrix elements for the states within the four- 
spinon continuum, at q = 0. We can simplify the problem 
by means of a very crude approximation and obtain a 
result which, although not rigorous, yet is meaningful for 
a comparison with the experiment. Let us replace in Eq. 
P2| the quantity (_E„ — Eg) with the average energy value 
of the four spinon continuum: Es — 7tJi\. In this way, as 
in the case of the CT excitations, we can write: 



^a,{u;)<h; = ^^^EsgsiT), 



(27) 



where gs{T) is the spin-correlation function defined as: 
gsiT) ^ A{{hl) ~ {hs)^)T . (28) 

One has to note that a non zero value of gs re- 
quires that the total Hamiltonian Ht and the interaction 
Hamiltonian Hs do not commute with each other. This 
condition would be realized including, in Ht or in ^,5, 
nnn interaction terms which would, unfortunately, com- 
plicate the problem considerably. However, a finite value 
of gs is obtained also in case of an AF broken symmetry 
of the ground state. We can then estimate an upper limit 
for gs starting from a long range Neel ordered state, i.e. 
\g) — Y[j \L2j,'\L2j+i.i)- Over this ground state we can 
calculate gs ^ N. On the other hand, for a random 
orientation of spins we would obtain gs = 0. The real 
spin configuration of the system is of course something 



in between these two extreme cases and the value of gs 
depends on the details of the many-body wave function 
of the spins. In particular, it depends on the probability 
of having fragments of three neighboring spins ordered 
antiferromagnetically. Therefore, gs is higher the lower 
is the temperature and reaches for T — > its maximum 
value, which would be strictly smaller than gs ^ N in 
a truly ID system. In relation to the interpretation of 
the optical conductivity data of a'-NaV205 within the 
charged magnon model, it is important to note that, in a 
dimerized singlet (or triplet) configuration of the single 
two-leg ladder, we have gs = 0. 

In conclusion, the intensity of the spin fluctuations rel- 
ative to the CT excitations in terms of effective charges, 
for the single two-leg ladder system, is given by: 
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On the basis of the analysis presented above, we can state 
that the maximum limiting value for this quantity, in the 
limit [/— >oo and with Jn given by Eq. n^, is: 
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IV. OPTICAL SPECTROSCOPY 

We investigated the optical properties of a'-NaV205 
in the frequency range going from 30 to 32 000 cm~^. 
High-quality single crystals were grown by high tempera- 
ture solution growth from a vanadate mixture flux.Lj The 
crystals, with dimensions of 1, 3, and 0.3 mm along the a, 
6, and c axes, respectively, were aligned by conventional 
Laue diffraction, and mounted in a liquid He flow cryo- 
stat to study the temperature dependence of the optical 
properties between 4 and 300 K. Reflectivity measure- 
ments, in near normal incidence conflguration (^^^10°), 
were performed on two different Fourier transform spec- 
trometers: A Bruker IFS 113v, in the frequency range 
going 20-7000 cm^^, and a Bomem DA3, between 6000 
and 32 000 cm^^. Polarized light was used, in order to 
probe the optical response of the crystals along the a and 
b axes. The absolute reflectivity was obtained by cali- 
brating the data acquired on the samples against a gold 
mirror, from low frequencies up to 15 000 cm^^, and an 
aluminum mirror for higher frequencies. The optical con- 
ductivity was calculated from the reflectivity data using 
Kramers-Kronig relations. 

Let us now, as an introduction to what we will dis- 
cuss in detail in the following sections, describe briefly 
the main features of the optical spectra of a'-NaV205, 
over the entire frequency range we covered with our ex- 
perimental systems. In Fig. |g we present reflectivity and 
conductivity data of a'-NaV205 at 300 K for E\\a (i.e.. 
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FIG. 5. Optical spectra of a'-NaV205 at 300 K for E \\ a 
(i.e., 1. to the chain direction) and E \\ b (i.e., || the chain 
direction), in the frequency range going from 30 to 17 000 
cm~^. Panel (a), and (b) show the results for reflectivity, and 
optical conductivity, respectively. In the inset of panel (b) an 
enlarged view of ai{uj) from 40 to 3000 cm~^ is presented. 



_L to the chain direction) and E\\b (i.e., || the chain di- 
rection), in the frequency range going from 30 to 17000 
cm~^. We can observe a very rich and peculiar exci- 
tation spectrum, characterized by strong anisotropy for 
direction parallel and perpendicular to the V-0 chains. 

For E II chain, the far-infrared region {oj < 800 cm~^) is 
characterized by strong optical phonon modes, with no 
background conductivity, as expected for an ionic insu- 
lator. At higher frequencies a strong electronic absorp- 
tion (peaking at '--^ 9000 cm~^) is directly observable in 
both reflectivity and conductivity. Moreover, in the mid- 
infrared region (between 800 and 5000 cm~^) we detected 
a very weak continuum of excitations, as shown by the 
finite value of the conductivity (see inset of Fig. |^) and 
by the refiectivity which is continuously raising in that 
frequency range (Fig. 0a). 

On the other hand, for i?_L chain, a weak broad band of 
optical absorption, additional to the phonon modes, was 
detected in the far-infrared region. This is directly recog- 
nizable in the reflectivity spectrum (Fig. Ha): Reducing 
the frequency of the incoming light, contrary to what ob- 
served for E II chain, we can observe an overall increase in 
reflectivity which masks the contribution due to the sharp 
phonon lines. At the same time, the interference fringes. 



due to Fabry-Perot resonances, which were dominating 
the optical reflectivity below 140 cm"-'^ for E\\ chain, are 
almost completely suppressed for E _L chain, indicating 
a stronger absorption for the latter polarization. In the 
optical conductivity, a continuum, with a broad maxi- 
mum in the far infrared (see inset of Fig. ^d), extends 
from very low frequencies up to the electronic excitation 
at ~ 7500 cm~^ ^ 1 eV, which is now much more intense 
than the similar one detected for £■ || chain. Moreover, 
along the direction perpendicular to the V-0 chains, we 
can observe an additional shoulder in the conductivity 
spectrum, between 3000 and 5000 cm~^. 

Finally, for both polarizations, an absorption edge due 
to the onset of charge-transfer transitions was detected 
at ^ 3 eV (not shown) , in agreement with Ref. |35| . 



A. Phonon Spectrum: High Temperature Phase 

We will now concentrate on the phonon spectra of 
the high temperature undistorted phase of a'-NaV2 05. 
We will try to assess the symmetry issue by compar- 
ing the number and the symmetry of the experimen- 
tally observed optical vibrations to the results obtained 
from the group theoretical analysis for the two differ- 
ent space groups proposed, for a'-NaV2 05, in the undis- 
torted phase. In Fig. o we present the reflectivity data 
for E perpendicular and parallel to the V-0 chains (i.e., 
along the a and b axes, respectively), up to 1000 cm~^, 
which covers the full phonon spectrum for these two crys- 
tal axes. In this respect, different is the case of the c 
axis: Along this direction a phonon mode w as de tected 
at ~ 1000 cm^^, as we will show in section IV C while 
discussing the results reported in Fig. O. In Fig. o the 
data are shown for T= 100 K, temperature which is low 
enough for the phonons to be sufficiently sharp and there- 
fore more easily detectable. At the same time, it is far 
from Tc, so that we can identify the phonon spectrum of 
the undistorted phase without having any contribution 
from the one of the low temperature distorted phase. 

In Fig. ra three pronounced phonon modes are observ- 
able for £'11 chain (wto « 177, 371, and 586 cm-\ at 100 
K), and three for E _L chain (cjto ~ 137, 256, and 518 
cm~^, at 100 K). It has to be mentioned that the fea- 
ture at 213 cm^^ along the a axis, which looks like an 
antiresonance, is a leakage of a c-axis modcEj due to the 
finite angle of incidence of the radiation on the sample 
(^'--^lO"), and to the use of p-polarization (therefore with 
a finite £^11 c component) to probe the a-axis optical re- 
sponse. Analyzing in detail the spectra few more lines 
are observable, as shown in the insets of Fig. 0. Along 
the b axis a fourth phonon is present at 225 cm~^: Be- 
ing the sample a 300 fiui thick platelet characterized by 
some transparency in this frequency region, Fabry-Perot 
interferences with a sinusoidal pattern are measurable. 
At 225 cm^"'^ the pattern is not regular anymore, indi- 
cating the presence of an absorption process. From the 
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FIG. 6. Reflectivity of a'-NaV205 in the undistorted phase 
at 100 K. The spectra are shown for E perpendicular and par- 
allel to the chain direction in panel (a) and (b), respectively. 
An enlarged view of the frequency regions, where very weak 
optical phonons were detected, is given in the insets. 



temperature dependence, it can be assigned to a lattice 
vibration. Similarly, three more phonons are detected 
along the a axis at 90, 740 and 938 cm~^. 

In conclusion, 6 a- axis and 4 &-axis phonons were de- 
tected on a'-NaV2 05, in the high temperature undis- 
torted phase. The resonant frequencies of these vibra- 
tional modes are summarized, for two different temper- 
atures, in table |. These results compare better with 
the group theoretical analysis for the centrosymmetric 



space group Pmmn (see Eq. ^ 

brational modes for E||a and E' 

the analysis for the noncentrosymmetric P2imn. In fact, 

in the latter case 15 and 7 phonons ought to be expected 



, which gave 7 and 4 vi- 
b, respectively, than with 



Pol. 


T(K) 


ujTO (cm ^) 


E\\a 


4 
100 


91.12 139.15 255.49 517.30 740.62 938.21 
89.67 136.81 255.82 517.74 740.49 938.49 


E\\b 


4 
100 


177.46 224.67 373.91 587.20 
177.18 225.20 371.02 585.88 



TABLE I. Resonant frequencies of the optical vibrational 
modes characteristic of the undistorted phase of a'-NaV205, 
obtained by fitting the optical phonons with Fano line shapes. 
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FIG. 7. Temperature dependence of resonant frequency 
shift, damping, and oscillator strength (obtained by fitting 
the data with Fano line shapes), for the most intense optical 
phonons of the undistorted phase of a'-NaV205. 



(see Eq. ||). Therefore, is. 
structure redeterminationJlj" 



jreement with the x-ray 
1 also the optical investi- 
gation of the phonon spectra indicates the centrosym- 
metric space group Pmmn as the most probable one for 
a'-NaV205, at temperatures above Tc- 

However, P2imn is in principle not completely ruled 
out as a possible space group on the basis of the optical 
phonon spectra. In fact, group theory gives information 
only about the number of phonons that one could ex- 
pect but not about their resonant frequencies and, in 
particular, their effective strength which, in the end, de- 
termines whether a mode is detectable or not. The only 
really conclusive answer would have been to detect more 
modes than allowed by symmetry for the space group 
Pm,m,n. In the present situation, it may be possible that 
the correct space group is still the non-centrosymmetric 
P2imn and that some of the phonons have escaped de- 
tection because the frequency is too low for our setup or 
due to a vanishing oscillator strength. In this respect, 
we have to stress that, indeed, the oscillator strength of 
the additional modes expected for the noncentrosymmet- 
ric P2imn would have to be small Jaecause, as indicated 
from the x-ray diffraction analysis ,Ej if a deviation from 
centrosymmetry is present, it is smaller than 0.03 A. We 
have also been performing a lattice dynamical calcula- 
tion for both space groups which shows that, removing 



the center of inversion and changing the valence of the V 
ions (i.e., from the uniform V^'^"*" for all sites to an equal 
number of V^^ and V^+ sites), the number of modes 
characterized by a finite value of the oscillator strength 
does not increase, even if a strictly long range charge or- 
dering in ID chains of V^+ and of V^+ ions is assumed. 
It has to be mentioned that many independent inves- 
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tigations of the Raman anc - 
a'-NaV205, were reported, 
not immediately a perfect agreement between all the dif- 
ferent data. As a matter of fajCt . in the first experimental 
studies of lattice vibrations,E30 the optical spectra there 
presented were interpreted as an evidence for the non- 
centrosymmetric space group P2imn. At this stage, we 
believe that the results of our investigation performed on 
high quality single crystals, give a very complete picture 
of the infrared vibrational modes of a'-NaV205, in the 
undistorted phase, for E\\a and E\\b. Moreover, our re- 
sults are completely confirmed by those recently obtained 
by Smirnov et aZHjH (who could measure also the c-axis 
phonon spectrum) , on samples of different origin. 

Before moving on to the phonon spectra characteristic 
of the low temperature distorted phase of a'-NaV205, we 
would like to discuss the temperature dependence, from 
4 to 300 K, of the lattice vibrations which are already 
present in the uniform phase. In order to extract, from 
the experimental data, information about their parame- 
ters, we had tait the data using Fano line shapes for the 
phonon peaks a In fact, most of them, in particular along 
the a axis, are characterized by a strong asymmetry indi- 
cating an interaction with the underlaying low frequency 
continuum. In this case the symmetrical Lorcntz line 
shape was not suitable and the asymmetrical Fano profile 
had to be used: The latter model contains an additional 
parameter which takes care of the degree of asymmetry 
of an excitation line, due to the coupling between that 
discrete state and an underlaying continuum of states. 

The results for the percentage change in resonant fre- 
quency, the damping, and the oscillator strength of the 
most intense optical phonons of the undistorted phase of 
a'-NaV205, are plotted in Fig. 0, between 4 and 300 K. 
One can observe that some of the modes show a sudden 
change in their parameters between 40 and 32 K, i.e., 
upon reducing the temperature below Tc = 34 K. In par- 
ticular, the oscillator strength of the phonons at 371 and 
586 cm^^, along the b axis, decreases across the phase 
transition, suggesting either a transfer of spectral weight 



Pol. (LT) 


LOTO (cm ^) 


Ella 


101 127 147 199 — 362 374 410 450 717 960 
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102 128 — — 234 362 — 410 — 718 960 



TABLE II. Resonant frequencies, at 4 K, of the zone 
boundary folded modes detected on a'-NaV2 05 in the low 
temperature (LT) phase, with i?_L chain and E \\ chain. 
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FIG. 8. Reflectivity spectra of Qf'-NaV205 measured for 
_E_L chain (a) and E ||chain (b) below (4 K) and just above (40 
K) the phase transition. Along both axes, new phonon lines 
activated by the phase transition are observable for T <Tc. 



to zone boundary folded modes activated by the phase 
transition or to electronic degrees of freedom. We will 
see later, in the course of the paper, that for the chain 
direction the reduction of the phonon oscillator strength 
will be compensated mainly by the strength gained by 
the zone boundary folded modes activated by the phase 
transition. On the other hand, for i?_L chain, the change 
in spectral weight of the electronic excitations will appear 
to be particularly important. 



B. Phonon Spectrum: Low Temperature Phase 

Upon cooling the sample below Tc=34 K, significant 
changes occur in the optical phonon spectra as one can 
see in Fig. g|, where we present the reflectivity spectra 
of a'-NaV205 measured, for i5_L chain (a) and -E|| chain 
(b), below (4 K) and just above (40 K) the phase tran- 
sition. Contrary to the case of CuGeOa, where we could 
clearly observe in ceflectivity only a single, very weak 
additional phonon,c3L3 10 new lines are detected for 
£;_L chain (wto~101, 127, 147, 199, 362, 374, 410, 450, 
717, and 960 cm^^), and 7 for £;|| chain (wto ~ 102, 128, 
234, 362, 410, 718, and 960 cm^^), some of them showing 
a quite large oscillator strength. The stronger intensity 
of the lines, with respect to the case of CuGeOs, suggests 
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FIG. 9. Detailed temperature dependence of some of the 
zone boundary folded modes observed in the reflectivity spec- 
tra, for T<Tc, with _E± chain (a,b) and E\\ chain (c,d). 



larger lattice distortions in a'-NaV205, on going through 
the phase transition or, alternatively, a larger electronic 
polarizability. The resonant frequencies (at T=4 K) of 
all the detected zone boundary folded modes are sum- 
marized in table O. The observation of many identical 
or almost identical frequencies for both axes, which is 
an important information for a full understanding of the 
structural distortion, is an intrinsic property of the low 
temperature phase. Experimental errors, like polariza- 
tion leakage or sample misalignment, are excluded from 
the well-defined anisotropy of the phonon spectrum of 
the undistorted phase. A more detailed discussion of this 



point will be presented in section [V C , in relation to the 
reflectivity data reported in Fig. nU 

An enlarged view of few of the frequency regions where 
zone boundary folded modes were detected, is given in 
Fig. ^ We observe the gradual growing of the these 
modes, whose oscillator strength increases the more the 
temperature is reduced below Tc=34 K. Moreover, the 
reflectivity data acquired for _E_L chain [panel (a) in Fig. 
II and panels (a) and (b) in Fig. g] show that also the un- 
derlaying continuum has a considerable temperature de- 
pendence: The reflectivity is decreasing, over the entire 
frequency range, upon reducing the temperature from 
40 to 4 K. In particular, in Fig. 0a, we see that very 
pronounced interference fringes, due to Fabry-Perot res- 
onances, are present at 4 K for uj < 136 cm~^, indicating 
a particularly strong reduction of the absorption in that 
frequency region. This effect is, in our opinion, related to 
the opening of the magnetic gap in the excitat ion sp ec- 
trum, and will be discussed in detail in section [VD. 



Of all the detected folded modes, for experimental rea- 
sons (i.e., higher sensitivity of the detector and therefore 
higher accuracy of the results in this frequency range) , we 
carefully investigated the ones at 718 (E \\ a, b) and 960 
cm~^ (E\\a), for several temperatures between 4 and 40 
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FIG. 10. Temperature dependence of the change in reso- 
nant frequency (a) and damping (c), of the normalized oscil- 
lator strength (d), and of the Fano asymmetry parameter (b), 
for the folded phonons observed below the phase transition at 
717 and 960 cm~^ with Ella, and at 718 cm"^ with £116. 



K. We fitted these phonons to a Fano profile because the 
718 cm~^ modes show an asymmetrical line-shape. The 
results of the percentage change in resonant frequency 
and damping, of the normalized oscillator strength, and 
of the Fano asymmetry parameter are plotted versus tem- 
perature in Fig. n3. We use here an asymmetry param- 
eter Q (rad) defined in such a way that the larger the 
value of 8, the stronger is the asymmetoz^of the line: A 
Lorentz line-shape is recovered for 8=0.CJ 

Let us now start by commenting on the results ob- 
tained for the oscillator strength: In Fig. [lC| d we see that 
5 has a similar behavior for the three different lines. How- 
ever, the 960 cm"^ peak vanishes at Tc whereas the two 
718 cm^^ modes have still a finite intensity at r=40 K 
and, as a matter of fact, disappear only for T> 60-70 K. 
At the same time, the line-shape of the 960 cni^-'^ mode is 
perfectly lorentzian at all temperatures, whereas the two 
other phonons show a consistently increasing asymmetry 
for T > 32 K. Also the increase in resonant frequency 
(Fig. [10|a) and the reduction of the damping (Fig. 10:) 
are much more pronounced for the two 718 cm~^ modes. 
From these results we conclude that the second-order 
character of the phase transition is nicely shown by the 
behavior of 5 for the 960 cm^^ folded mode. On the other 
hand, pretransitional fluctuations manifest themselves in 
the finite intensity of the 718 cm"'^ modes above Tc, and 
in the extremely large value of and 7 for these lines 
(7'--^35 cm~^ at 40 K). In other words, the lattice distor- 
tion is already taking place in the system, at temperature 
much higher than Tc, but with a short range character, 
whereas a coherent long range distortion is realized only 
below the phase transition temperature. In this respect, 
it is worth mentioning that similar pretransitional fiuctu- 
ations, below 70 K, have been observed also in the course 
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FIG. 11. Reflectivity of a'-NaV205, at 6 K, for several 
orientations of the electric field of the light in the a-b plane (as 
sketched in panel b, with A(f> = 10°). The values = 0°, 180° 
correspond to _E± chain, and = 90° to E\\ chain. 



of a study of the propagation of uLtrasonic waves along 
the chain direction of Q:'-NaV205,Ej ancL-yery recently, 
in x-ray diffuse scattering measurements. Ell 



C. Symmetry of the Lattice Distortion 

We saw in the previous section that many of the folded 
zone boundary phonons, activated by the phase transi- 
tion, show the same resonant frequency for E\\a and E\\b 
(see table ||), even though they are characterized by a 
different oscillator strength along the two different axes 
(e.g., S'~0.021 for E\\a and 5-0.014 for E\\b, for the 718 
cm~^ modes at 4 K). This is quite a surprising result be- 
cause the phonon spectrum of the undistorted phase has 
a well defined anisotropy, at any temperature, with dif- 
ferent resonant frequencies for vibrations polarized along 
the a and b axes (see table |) . Similar anisotropy should 
then be present below the phase transition if the system 
is still orthorhombic, even for the folded modes. 

To further check these findings, we performed reflec- 
tivity measurements on a single crystal of a'-NaV205, at 
the fixed temperature of 6 K, rotating t he p olarization of 
the incident light of an angle (j) (see Fig. |ll|). The sample 
was aligned in such a way to have the electric field of the 



light parallel to the a axis of the crystal, for = 0° and 
180°, and parallel to the b axis, for = 90°. For these 
polarizations of the light the obtained results are, obvi- 
ously, identical to those presented in Fig. ||. However 
in Fig. O the data are shown in a different frequency 
range, extending up to 1050 cm^^. We can then ob- 
serve at — 1000 cm^^ along the a axis, which again was 
probed with p-polarizcd light j— the antiresonance due to 
the leakage of a c-axis phonon. e£| Starting from = 0° and 
turning the polarization in steps A0= 10°, we observe a 
gradual decrease of the a axis contribution and, at the 
same time, an increase of the b axis contribution to the 
total reflectivity. Moreover, the set of curves plotted in 
Fig. O satisfy the following relations: 

i?(45°) + i?(-45°) = i?(0°) + i?(90°) , (31) 

i?(0) = i?(O°)cos2(0) + i?(9O°)sm2(0) . (32) 

This suggests that the system, even in the distorted 
phase, has the optical axes along the crystal directions a, 
6, and c (on the contrary, for a triclinic or a monoclinic 
structure there would be no natural choice of the opti- 
cal axes and Eq. ^ would not be satisfied). However, 
the low temperature folded mode at 718 cm~^ can be 
observed for all possible orientations of the electric field 
of the incident radiation (see Fig. W^. Because identical 
resonance frequencies, along the a and the b axes of the 
crystal, were observed for six of the folded zone boundary 
modes, we cannot believe in an accidental coincidence. 

A possible explanation for the observed effect is as 
follows. In the low temperature phase a non-magnetic 
ground state is realized, and true singlets are present in 
a'-NaV205. The most probable situation is that the sin- 
glets are formed, within the ladders, over the plaquettes 
defined by four V ions. In this case, as shown in Fig. O, 
within one ladder there is an alternation of 'empty' pla- 
quettes and plaquettes containing a singlet (represented 
by a grey ellipsoid in Fig. |l^). Let us now assume that 
in adjacent ladders the charge is arranged in an oblique 
charge ordering pattern (OCOP), as shown in the left- 
hand side of Fig. |l^. The OCOP is consistent with the 
doubling of the lattice constant observed, along the u 
axis, in low temperature x-ray scattering experimentsB 
Moreover, two different kinds of V ions are identifiable 
which could possibly explain the-results obtained in low 
temperature NMR expcrimcntstj (for a more detailed 
discussion see section 0). In fact, there are V ions adja- 
cent to empty plaquettes (white circles) and V ions ad- 
jacent to plaquettes containing a singlet (black circles), 
arranged in zig-zag patterns on each ladder. Note that 
in our picture the non-equivalence of plaquettes along 
the 6-axis, indicated by the gray ellipsoids in Fig. |lj, 
is really a chicken-and-egg problem: Either the singlet 
formation occurs first, causing an alternation along b of 
non-equivalent plaquettes, which, in combination with 
the OCOP along a gives rise to a zig-zag charge order- 
ing. Alternatively, the zig-zag charge order occurs first. 
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FIG. 12. Possible charge ordering pattern in a'-NaV205, 
below Tc- The grey ellipsoids represent the electronic charge 
involved in the formation of singlets, within the ladders, over 
every second plaquette formed by four V ions. On adjacent 
ladders, the charge is assumed to define an oblique charge or- 
dering pattern (OCOP; see, e.g., left-hand side). Two differ- 
ent kinds of V ions are identifiable: V ions adjacent to empty 
plaquettes (white circles) and V ions adjacent to plaquettes 
containing a singlet (black circles), arranged in zig-zag pat- 
terns on each ladder. If domains with different orientation 
for the OCOP are present in the system (left and right-hand 
sides), the zig-zag ordering of the different kind of V ions is 
destroyed in correspondence of a domain wall (dashed line) . 



which, in combination with OCOP causes the alterna- 
tion of non-equivalent plaquettes along b. The resulting 
charge modulations are only small deviations from aver- 
age occupation, which we represented graphically by the 
different grey shades. 

At this stage, the low temperature phonon spectra 
should still be characterized by a well defined anisotropy, 
with distinct eigenfrequencies for different polarizations 
of the electric field of the incoming radiation. However, 
there is no reason for the polarization of the folded zone 
boundary modes to be precisely along the a or the b axes. 
In fact, they are polarized with a finite angle with respect 
to the crystal axes (e.g., along directions parallel and per- 
pendicular to the diagonal long-range singlet-pattern) . If 
domains are present in the system, with opposite orienta- 
tion for the OCOP (left and right hand sides of Fig. n2), 
each domain will still be characterized by a phonon spec- 
trum with a well defined polarization, as just discussed. 
However, in presence of such domains the phonon spectra 
measured in the experiments correspond to an average 
over many domains with the two possible orientations 
randomly distributed. As a result, almost no anisotropy 
would be found for light polarized along a and b axes. 

A short wavelength alternation of descending and as- 
cending OCOP domains (also a. 2a x 2b supercell) cor- 



responds closeli' to the structure recently reported by 
Liidecke et alJtB A similar structure was also considered 
by Riera and Poilblanc as one of the many possible real- 
izations of charge density wave in this systemO 



D. Optical Conductivity 

Very interesting information and a deeper understand- 
ing of the symmetry of Q;'-NaV205 are obtained from 
the detailed analysis of the optical conductivity data, in 
particular when the electronic excitations are considered 
and compared to the result relative to CuGeOs (see Fig. 
O). On the latter comnaund we observed sharp phonon 
lines below 1000 cm~^,cErc3 multiphonon absorption at 
~1500 cm~^, very weak [note the low values of ai{uj) 
in the inset of Fig. [l^a] jjhonon- assisted Cu d-d tran- 
sitions at ~ 14000 cm~^,E3 and the onset of the Cu-0 
charge-transfer (CT) excitations at ~27 000 cm~^. On 
a'-NaV205, besides the phonon lines in the far-infrared 
region, we detected features that are completely absent in 
CuGeOs: A strong absorption peak at ~8000 cm~^ and, 
in particular for E _L chain, a low-frequency continuum 
of excitations (see inset of Fig. [l3|b). Our goal is to in- 
terpret the complete excitation spectrum for a'-NaV205, 
trying, in this way, to learn more about the system on a 
microscopic level. In order to do that, we have first to 
identify possible candidates for the excitation processes 
observed in the experimental data. 

On the basis of intensity considerations, the peak at 
~8000 cm^^^leV, in the optical conductivity of a'- 
NaV2 05, has to be ascribed to an optically allowed exci- 
tation. Contrary to the case of the absorption detected 
at ~ 14 000 cm^^ on CuGeOs, it cannot be interpreted 
just as a d-dexciton on the transition metal ion (i.e., V in 
the case of Q;'-NaV205). In fact, this kind of transition is 
in principle optically forbidden and only weakly allowed 
in the presence of a strong parity-broken crystal field or 
electron-phonon coupling. We have now to recall one of 
the main consequences of the structural analysis which 
suggested as a most probable space group, for a'-NaV205 
in the undistorted phase, the centrosymmetric Pmmn: 
The existence of only one kind of V per unit cell, with 
an average valence of -1-4.5. Taking into account that 
the hopping parameter along the rungs t±, is approxi- 
mately a factosjpttwo larger than the one along the legs of 
the ladder i|| ,L3e3 all the other hopping amplitudes being 
much smaller (including the nn t^y), a'-NaV205 can be 
regarded as a 1/4-filled two-leg ladder system, with one 
V d-electron per run|g-|Shared by two V ions in a molec- 
ular bonding orbital.E3 In this context, the most obvious 
candidate for the strong optical absorption observed at 
1 eV is the on-rung bonding-antibonding transition, with 
a characteristic energy given by 2tj_. Hossiever, on the ba- 
sis of LMTO band structure calculationE^ and exact idiag- 
onalization calculation applied to finite size clusters)^ t± 
was estimated to be '--^ 0.30-0.35 eV, which gives a value 
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FIG. 13. Optical conductivity, at 300 K, of CuGeOs (a) 
and a'-NaV2 05 (b), for E parallel and perpendicular to the 
Cu-O or to the V-O chains. Inset of panel a: enlarged view of 
o-i(u) of CuGeOs, from 800 to 30 000 cm"^Inset of panel (b): 
enlarged view of ai(Lo) of a'-NaV2 05, from 40 to 3000 cm~^. 
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FIG. 14. Far-infrared optical conductivity of Q'-NaV205, 
for E _L chain (a) and E || chain (b). Along both axes new 
phonon lines, activated by the phase transition, are present 
for T < Tc = 34 K. Moreover, for E _L chain, the temperature 
dependence of the low frequency continuum is observable. 



for the bonding-antibonding transition lower than the ex- 
perimentally observed 1 eV. This energy mismatch, and 
the strong dependence of the intensity of the 1 eV peak 
on the polarization of the light, are important pie ces of 
information and will be discussed further in section IV E| . 

Let us now turn to the infrared continuum detected for 
i?_L chain (see inset of Fig. ^SJd). We already showed , 
discussing the low temperature reflectivity data (see Fig. 
pi and P), that this broad band of optical absorption has 
a considerable temperature dependence: The overall re- 
flectivity was decreasing upon reducing the temperature 
and, for T <Tc, Fabry-Perot resonances appear at fre- 
quencies lower than 136 cm~^, indicating a particularly 
strong reduction of the optical absorption in that fre- 
quency region. These effects become even more evident 
when we consider the low temperature optical conduc- 
tivity calculated from the reflectivity data via Kramers- 
Kronig transformations (see Fig. O). However, we have 
to keep in mind that performing this kind of calculation, 
we will obtain unphysical oscillations for the optical con- 
ductivity in the frequency regions where reflectivity is 
dominated by interference fringes. 

In Fig. [l4p we can observe, for E \\ chain, only the 
sharpening of the phonon lines and the appearance of 
the zone boundary folded modes. On the other hand. 



for E _L chain we can observe, in addition to the folded 
modes below Tc, first the increase of the continuum in- 
tensity upon cooling the sample from 300 to 40 K and, 
subsequently, a reduction of it from 40 to 4 K. In the 
frequency region below the 139 cm""'^ phonon, whereas a 
reduction is observable in going from 300 to 40 K, not 
much can be said about the absolute conductivity of the 
4 K data, because of the interference fringes. 

In order to evaluate the low temperature values of 
conductivity more precisely, we can fit the interference 
fringes directly in reflectivity. In fact, fitting the period 
and the amplitude of them, we have enough parameters 
to calculate right away the real and imaginary part of 
the index of refraction and, eventually, the dynamical 
conductivity cti . The so obtained values are plotted ver- 
sus the temperature in Fig. [l5|a. The data for 100, 200, 
and 300 K were obtained via Kramers-Kronig transfor- 
mations, whereas those for 4 and 32 K from the fitting of 
the interferences in the range 105-125 cm~^. The 40 K 
point, as at that temperature very shallow fringes start to 
be observable, was calculated with both methods, which 
show a rather good agreement in this case. As a results 
we can observe (see Fig. 15 1) a pronounced reduction 
of cTi, across the phase transition, corresponding to the 
opening of a gap in the optical conductivity. 
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FIG. 15. Panel (a): Temperature dependence of the optical 
conductivity at lo = 115 crn"^, for i?_L chain. The values are 
obtained from Kramers-Kronig analysis, for T > Tc, and by 
fitting the interference fringes in reflectivity, for T <Tc. The 
inset shows an enlarged view of cri(aj), from 110 to 190 cm~^. 
Panel (b): Oscillator strength of the low frequency electronic 
continuum, for _E± chain, plotted versus temperature. 



of the underlying continuum has decreased, a sharp and 
rather symmetrical line shape is recovered. 

An additional information we can extract from the con- 
ductivity spectra is the value of the integrated intensity, 
and its temperature dependence, for the infrared broad 
band of optical absorption detected for E ± chain. The 
oscillator strength, obtained by integrating (Ti{lu) (with 
phonons subtracted) up to ^800 cm~^, is displayed in 
Fig. n5p. It increases upon cooling down the sample 
from room temperature to Tc, and rapidly decreases for 
T < Tc- The possible meaning of this behavior will be 
further discussed in the next section. 

As a last remark, it is important to mention that low 
temperature measurements of the dielectric constant, in 
the microwave frequency range, wei re , r ecently performed 
on a'-NaV205 by different groups.E^EZl The reported re- 
sults show no anomaly, across the phase transition, for 
the real part of the dielectric constant along the b axis. 
On the other hand, along the a axis, a pronounced de- 
crease below Tc was found. The observed behavior, as 
a matter of fact, is in very good agreement with the 
temperature dependence of the oscillator strength of the 
low frequency continuum present, in the optical spectra, 
along the a axis (Fig. [l5|b). This comparison is meaning- 
ful because the contribution of all possible optical exci- 
tations to the static dielectric constant is just given by 
their oscillator strengths. The very good agreement of 
the different results proves, as already suggested by the 
authors of Ref. 56 and p% that the change in the very 
low frequency dielectric constant is related to electronic 
(magnetic) degrees of freedom. 



It is difficult to evaluate precisely the gap size because 
of the phonon at 139 cm~^ (inset of Fig. [l5|a). How- 
ever, we can estimate a lower limit for the gap of ~ 136 
cm~^ (17±3 meV). The important information is that the 
value of the gap in the optical conductivity is approxi- 
mately twice the spin gap value observed in inelastic neu- 
tron scattering exjperimentaa and magnetic susceptibil- 
ity measurements, E3 for T< 34 K. This finding indicates 
that the infrared broad band of optical absorption has to 
be related to electronic degrees of freedom: Because the 
range of frequency coincides with the low-energy-scale 
spin excitations, the most likely candidates for this con- 
tinuum are excitations involving two spin flips. 

In relation to what we just discussed, it is interesting 
to note the effect of the opening of a gap, in the optical 
conductivity, on the phonon located at 139 cm"-'^. This 
can be clearly identified in the inset of Fig. HSa, where 
an enlarged view of (Ti{uj) at low frequency is given, for 
100 and 4 K. At 100 K this mode is characterized by a 
very pronounced anomaly of its line shape. As a matter 
of fact, fitting this phonon with a Fano line shape we 
obtained for the asymmetry parameter the value 6 = 
— 1.7 rad, indicating a strong coupling to the continuum. 
However at 4 K, once the gap is open and the intensity 



E. Charged Bi-Magnons 

In the previous section we identified possible can- 
didates for the 1 eV excitation (on-rung bonding- 
antibonding transition), and for the low frequency con- 
tinuum detected on Q;'-NaV205 for E _L chain (double 
spin flips). However, we also stressed that, for a sym- 
metrical 1/4 -filled two leg ladder system, the bonding- 
antibonding transition would be at 2f^'^0.7eV and not 
at the experimentally observed value of 1 eV. Moreover, 
also intensity and polarization (_E_L chain) of the contin- 
uum cannot be understood assuming the complete equiv- 
alence of the V sites required by the space group Pmmn. 
In fact, because of the high symmetry, no dipole mo- 
ment perpendicular to the V-0 chains would be associ- 
ated with a double spin flip process. Therefore, no opti- 
cal activity involving magnetic degrees of freedom would 
be detectable. Not even phonon-assisted-agp. excitations 
considered by Lorenzana and SawatzkytSEZfO (where the 
involved phonon is effectively lowering the symmetry of 
the system) could be a possible explanation for the ob- 
served spectra: The energy is too low for this kind of pro- 
cesses which, on the other hand, are probably responsi- 
ble for the weak mid- infrared continuum in the E\\ chain 
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spectra (see Fig. |l3| b and, in particular, its inset). 

A possible solution is to assume a charge dispropor- 
tionation on each rung of the ladders without, however, 
any particular long range charge ordering (see Fig. |l6| ). 
The latter is, in fact, most probably excluded by the 
x-ray diffraction results and by the phonon spectra of 
the undistorted phase. It is worth mentioning that the 
assumption of charge disproportionation has been con- 
firmed by calculations of the optical conductivity hu 
exact-diagonalization technique on finite size clustersE3 
Only in this way it was possible to reproduce the energy 
position of the 1 eV peak, and the anisotropy between a 
and b axes (see Fig. |3|). We can now try to interpret 
the optical data on a'-NaV205 in the framework of the 
charged- magnon model developed in section III . 

Let us start from the analysis of the absorption band at 
1 eV. We can model the jth rung, formed by two V ions, 
with the Hamiltonian Hjq + Hj±, as in Eq. and 0, where 
the potential energy difference between the two sites A 
has the role of a charge disproportionation parameter: 
For the symmetric ladder A = 0. The splitting between 
the two bonding and antibonding eigenstates \L) and \R), 
given by Ect = a/A^ -|- 4i5_, corresponds to the photon 
energy of the optical absorption. We see that, if A ^ 0, 
the experimental value of 1 eV can be reproduced. 

A second crucial piece of information is provided by 
the intensity of the absorption. We can calculate the in- 
tegrated intensity of the 1 eV peak, in the conductivity 
data, and then take advantage of its functional expres- 
sion derived in section 111 (Eq. y). This way we obtain 
from the spectra \t± \ « 0.3 eV, value which is in very good 
agreement with those obtained from LMTO band struc- 
ture calculation,E3 and exact diagonalization calculation 
applied to finite size clusters. Efl Combining this number 
with EcT — ^S:Y, we obtain A«0.8eV. The correspond- 
ing two eigenstates have 90% and 10% character on either 
side of the rung. Therefore, the valence of the two V-ions 
is 4.1 and 4.9, respectively, and the optical transition at 
1 eV is essentially a charge transfer (CT) excitation from 
the occupied V 3d state at one side of the rung to the 
empty 3d state at the opposite side of the same rung. 

Let us now turn to the infrared continuum detected for 
i^_L chain. The presence of two V states (|L) and \R)) 
per spin, along with the broken left-right parity of the 
ground state on each rung, gives rise to the fascinating 
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behavior of the spin flips described in detail in section III 



for the 1/4 filled single two-leg ladder, and in a pictorial 
way for a'-NaV205, in Fig. O: If in a small fragment of 
the ladder each spin resides in a |L) orbital, with some 
admixture of \R), the inclusion of the coupling between 
the rungs, via the Hamiltonian i^y (Eq. ||), has no effect 
when the spins are parallel, due to the Pauli-principle. 
If the spins form an S" = state, the ground state gains 
some kinetic energy along the legs by mixing in states 
where two electrons reside on the same rung. In the limit 
[/— >cxD, these states have one electron in the \L) and one 
in the \R) state on the same rung. As a result, there is a 
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FIG. 16. Pictorial description of the electric dipole moment 
associated with a single spin-flip in a'-NaV205, when a local 
breaking of symmetry in the charge distribution on each rung 
is assumed. If three neighboring antiparallel spins are laying 
on the same leg of the ladder (top) , the charge of the j th elec- 
tron is asymmetrically distributed not only over the j th rung 
but also on the nn rungs (with opposite asymmetry) , because 
of the virtual hopping iy . If we flip its spin (bottom), the j th 
electron is confined to the j th rung, resulting in a reduced 
charge density on the right leg of the ladder. Therefore, there 
is a net dipole displacement (along the rung direction) of the 
antiparallel spin configuration compared to the parallel one. 



net dipole displacement of the singlet state compared to 
the triplet state: The spin-flip excitations carry a finite 
electric dipole moment along the rung direction. 

From the expression for the exchange coupling con- 
stant calculated in the limit C/— >oo (Eq. |l^), and using 
the values of t± aad A calculated from the optical data 
and |i||| ~0.2eV,EJ we obtain J|| k^SO meV. It is signifi- 
cant that this value is comparable to those reported on 
the basis of magnetic susceptibility measurements, which 
range from 38 to 48 meV (see Ref . ^ and 0, respectively) . 
Moreover, for the effective charge of the bi-magnon exci- 
tation process (Eq. ^^, we obtain qm/qe = 0-07. 

We can also compare the experimental intensity of the 
spin-fluctuations relative to the CT excitations to the 
maximum limiting value calculated for this quantity (Eq. 
pOI) . For the parameters obtained from the conductiv- 
ity spectra, the maximum relative intensity is ~0.GG14. 
The experimental value is ~0.0008, at T=300 K, in good 
agreement with the numerical estimate. 

We can now try to interpret the temperature depen- 
dence of the spin-fluctuation continuum shown in Fig. 
n5p: The oscillator strength increases upon cooling the 
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sample from room temperature to Tg, and rapidly de- 
creases for T<Tc. Within the charged-magnon model, 
the increase marks an increase in short range AF corre- 
lations of the chains: In fact, the intensity of the spin- 
fluctuations is proportional to the spin-spin correlation 
function gs{T) defined in Eq. |2^. Below Tc, nn spin- 
singlet correlations become dominant, and gs is sup- 
pressed, along with the intensity of the spin-fluctuations. 
In conclusion, by a detailed analysis of the optical con- 
ductivity, we provided a direct evidence for a local break- 
ing of symmetry in the electronic charge distribution over 
the rungs of the ladders in a'-NaV205, in the high tem- 
perature undistorted phase. Moreover, we showed that a 
direct two-magnon optical absorption process is respon- 
sible for the low frequency continuum observed perpen- 
dicularly to the V-0 chains. 



V. DISCUSSION 

As we discussed in the course of this paper, the inter- 
pretation of the phase transition in a'-NaV205 is still 
controversial. On the one hand, the SP picture_was sug- 
gested by preliminary magnetic susceptibilityEJ and in- 
elastic neutron scattering measurements.cl On the other, 
x-ray diffraction redeterminations of the crystal struc- 
ture of the cpnippund, in the high temperature undis- 
torted phase,E2l l3 showed that the dimerization within 
well defined ID magnetic chains is not very probable in 
a'-NaV205. Moreover, the magnetic field dependence of 
Tc, expected for artaue SP system on the basis of Cross 
and Fisher theory|3'cl has not been observed.l3l3 

Additional— |information came recently from NMR 
experiments :ll3 Above Tc, only one kind of V site (mag- 
netic) was clearly detected, whereas no signature of non- 
magnetic V^^ was found. However, another set of mag- 
netic V sites which areJnvisible in the NMR spectrum 
could not be ruled outJij On the other hand, below the 
phase transition, two inequivalent sets of V sites were de- 
tected and assigned to V^"*" and V^"'' . These results were 
then interpreted as an indication of the charge-ordering 
nature of the phase transition, in Q;'-NaV205.E3 Subse- 
quently, many |tljeoretical studies have been addressing 
this possibility.E3il3 The main point to be explained is 
why a charge ordering would automatically result in the 
opening of a spin gap (with the same transition temper- 
ature). Among all the models, the more interesting seem 
to be those were a 'zig-zag' charge ordering on the lad- 
ders is realized. In this case, the ordering would imme- 
diately result in a nonmagnetic ground state, due to the 
formation of spin_singlets on nn V sites belonging to two 
different laddersE3 or, altecnatively, to the modulation of 
J|| within a single ladder .Ii3 

In the discussion of our spectra of the optical conduc- 
tivity, we showed that the infrared continuum, observed 
perpendicularly to the V-0 chains direction, could be at- 
tributed to direct optical absorption of double spin-flip 



excitations. However, in order to have optical activity 
for these processes and, at the same time, to explain the 
energy position and the polarization dependence of the 
on-rung bonding-antibonding transition at 1 eV, we had 
to assume a charge disproportionation between the two 
V sites on a rung. This assumption seejmSptp contrast the 
results of the Xpiay diffraction analysialHli^ and of NMR 
measurements. E3 As far as the former is concerned, it is 
obvious that if the charge disproportionation does not 
have a long range character it would not be detectable 
in an x-ray diffraction experiment, but would still result 
in optical activity of charged bi-magnon excitatkps. Re- 
garding the disagreement with the NMR results,E£l we saw 
that the existence of another kind of V sites cpuld not be 
completely excluded in the NMR spectrum. Ej Further- 
more, in comparing the optical to the NMR data, there 
could also be a problem with the different time scales of 
the two experimental techniques (i.e., 'fast' and 'slow' for 
optics and NMR, respectively). Therefore, if there are 
charge fluctuations in the system, as one would expect 
for a charge ordering phase transition, the experimental 
results could be different. 

From the analysis of the optical conductivity, we can- 
not really say which is the more appropriate picture for 
the phase transition in a'-NaV205. We saw that, indeed, 
the electronic and/or magnetic degrees of freedom are 
playing a leading role. This is shown, e.g., by the strong 
temperature dependence of the infrared spin fluctuation 
spectrum which, in turn, determines the aflnmalous be- 
havior of the microwave dielectric constant .E3'E3 However, 
the reduction of the oscillator strength of the infrared 
continuum for T < Tc, which we discussed as a conse- 
quence of nn spin-singlet correlations becoming dominant 
below Tc, could also result from a zig-zag charge ordering 
on the ladders. In this case, as for the perfectly symmet- 
rical 1/4-fllled ladder, the optical activity of double spin 
flips would simply vanish. 

In our opinion, for a better understanding of the phase 
transition, the pretransitional fluctuations observejd,with 
optical spectroscopy, ultcasonic wave propagation]^ and 
x-ray diffuse scattering ,Eil have to be taken into account. 
They probably reffect the high degree of disorder present 
in the system above the phase transition, which is also 
shown by the local and 'disordered' charge disproportion- 
ation necessary to interpret the optical conductivity. In 
this context, a charge ordering transition, where more de- 
grees of freedom than just those of the spins are involved, 
seems plausible. It would also explain the experimentally 
observed entropy reduction, which appearecLto be larger 
than what expected for a ID AF HB chain.cl 

Also the degeneracy of the zone boundary folded 
modes, activated by the phase transition, could be an 
important piece of information. However, at this stage 
of the investigation of the electronic and magnetic prop- 
erties of Q;'-NaV205, very important insights could come 
from the x-ray diffraction determination of the low tem- 
perature structure. Probably, only at that point it will be 
possible to fit all the different results together and come 
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up 



with the final interpretation of the phase transition. 



VI. CONCLUSIONS 

In this paper we have been discussing the high and 
low temperature optical properties of a'-NaV205, in the 
energy range 4meV-4eV. We studied the symmetry of 
the material, in the high temperature undistorted phase, 
comparing the different crvstal—atructurcs proposed on 
the basis of x-ray diffractioroE3"E3 to our infrared phonon 
spectra. We found that the system is better described by 
the centrosymmetric space group P7n7nn,piri|-a|greement 
with the recent structure redetermination.ll3~t3 

On the other hand, by a detailed analysis of the elec- 
tronic excitations detected in the optical conductivity, 
we provided direct evidence for a charged disproportion- 
ated electronic ground-state, at least on a locale scale: A 
consistent interpretation of both structural and optical 
conductivity data requires an asymmetrical charge dis- 
tribution on each rung, without any long range order. 
We showed that, because of the locally broken symme- 
try, spin-flip excitations carry a finite electric dipole mo- 
ment, which is responsible for the detection of charged 
bi-magnons in the optical spectrum for E _L chain, i.e., 
direct two-magnon optical absorption processes. 

By analyzing the optically allowed phonons at various 
temperatures below and above the phase transition, we 
concluded that a second-order change to a larger unit 
cell takes place below 34 K, with a fluctuation regime 
extending over a very broad temperature range. 

The charged-magnon model we have been discussing 
in this paper has been developed explicitly to interpret 
the peculiar optical conductivity spectra of Q;'-NaV205. 
However, our opinion is that the importance of this model 
goes beyond the understanding of the excitation spec- 
trum of just this particular system. We expect this model 
(obviously in a more general and rigorous form), to be 
relevant to many of the strongly correlated electron sys- 
tems, where the interplay of spin and charge plays a cru- 
cial role in determining the low energy electrodynamics. 

At this point we would like to speculate that a specific 
case where charged magnon excitations could be observed 
is the stripe phase of copper oxides superconductora£3 
and of other strongly correlated oxides .Ell In fact, even 
though the magnetism is quasi 2D in these systems, in 
the stripe phase a symmetry-breaking quasi ID charge 
ordering takes place, which could result in optical ac- 
tivity for the charged magnons. For instance, in the 
optical conductivity spectra o£|Lai.67Sro.33Ni04 below 
the charge ordering transition,E3 it is possible to observe 
some residual conductivity inside the energy gap, whose 
nature has not been completely understood so far. With- 
out further experimental and theoretical investigations it 
is of course not possible to make a strong claim. However, 
the reported features are suggestive of charged magnon 
excitations. In fact, for the value of J obtained from the 



two-magnon Raman scattering on this materialj£j direct 
two-magnon optical absorption is expected in the same 
energy range (i.e., from zero to ~0.2 eV) where tke resid- 
ual structure has been experimentally observed. Ej 

Pushing the speculation to its limit, one could imagine 
that even the local breaking of symmetry introduced by 
a single impurity in an AF could result in optical activ- 
ity for pure spin-flip excitations. One particular system 
where this situation is probably realized is Zn-substituted 
YBa2Cu306, where Zn is introduced in the Cu02 planes. 
In the cuprate parent compounds direct bi-magnon op- 
tical absorption is not allowed due to the inversion sym- 
metry; only phonon-assisted magnetic excitations are op- 
tically active, whenever the involved phonon is lowering 
the symmetry of the system. In fact, phonon assisted bi- 
magnon excitations in the mid-infrared region were de= 
tected on both YBa2Cu306 and YBa2Cu2.85Zno.i5060 
But on the Zn substituted system it is possible to observe, 
in the optical conductivity spectra, various features with 
frequencies coincident with those of pure double spin-flip 
excitations and, therefore, most probably ascrihable to 
direct charged bi-magnon absorption processes.E3 
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